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BLOCK PROGRAMMING 


A Valuable Aid to Simulation of Lumped-Parameter Physical Systems 


The analog computer, or differential analyzer, isa 
powerful tool for obtaining dynamic solutions to differ- 
ential equations. Generally the equations to be sol- 
ved relate to a physical system of particular interest 
to the computer programmer, in which case the pro- 
grammer wishes to simulate his physical system by 
means of an electrical network whose defining equa- 
tions are analogous (hence the name analog computer) 
to the subject physical system. It is desirable then, 
to arrange the computer mechanizationto have a one- 
to-one correspondence betweenthe problem board and 
the physical system, so that the programmer ‘'sees”’ 
his system rather than an abstract electrical network. 


The following discussion will illustrate a method of 
programming in which the physical elements are 
treated as blocks rather than as sets of differential 
equations. Certain useful basic building blocks will 
be developed, and methods for combining the basic 
blocks to form any desired physical system (simula- 
tion) will be discussed. Simulations generated by the 
methods to be described will have the advantage that 
any or all physical parameters appear as single con- 
trols in the computer so that any parameter may be 
changed without affecting any other setting. In addi- 
tion, any or all variables appear as computer outputs 
and may be individually monitored. 


A lumped-parameter system, by definition, consists 
of a number of interconnected discrete elements whose 
individual or collective responses to impressed forces 
or stimuli are of analytical interest. The forces and 
responses are related by a mathematical operator. 
Block programming starts with the selection of basic 
definition and rules. 


MATHEMATICAL OPERATOR (0) 


The operator, or combination of operators, describes 
the functional relationship between quantities. The 
most common operators are the following: 

Summer 


Constant Multiplier (potentiometer, commonly 
referred to as "pot"’) 


Inverter - -1 Sign Changer 
1 
Integrator - -> (La Place Transform Notation) 
Differentiator - s (This operator is never used in 
an analog simulation if it can 
be avoided, as itusually can). 
Variable Multiplier 
Arbitrary Function - f(x) 


Analytical Function Trigonometric, Log, 


Hyperbolic, etc. 


Inverse Operator - (0)? S, + are inverses 


ELEMENT 


The element is the basic unit of the block representation 
of a system. It hasa pair ofterminals or nodes, with 
a value associated with eachnode. It alsohas a trans- 
ference quantity between nodes, or ''through" the ele- 
ment whichis functionally related by some operator to 
the node-pair value. The transference quantity may 
also be referred to as the branch transference or trans- 
mission. 


The elements required to depict the passive lumped- 
parameters of most physicalsystems are very few in 
number. Inparticular, five basic elements will suffice 
for many linear electrical, mechanical, and thermal 
systems. These arethe Summer, Constant Multiplier, 
Sign Changer, Integrator, and Differentiator. 


BLOCK 


A block is the computer mechanization related toa 
physicalelement. The input/output variables of acom- 
puter block are voltages andcurrents, but the operator 
relationship is the same as the simulated element. 


Table I gives examples of single elements and their 
corresponding blocks. Table IIgives afew commonly 
occurring two-element combinations. Combinations 
of more than two elements can be generated as shown 
in the following examples. Usefultables of more com- 
plicated blocks or transfer-function simulations appear 
in many publications. 


DRIVING SOURCE 


A driving source is a source of energy or power for a 
network of elements. 


SYSTEM STABILITY 


A stable system is one whose responses are bounded 
(finite amplitude limit) for any finite input. An un- 
stable system is one whose responses are not bounded. 


NETWORK 


A network is an interconnected set of elements and 
sources. Elements may be joined together at their 
nodes, providing the connected nodes have common 
dimensions and value, and provided the transference 
quantities have common dimensions. 


This discussion concerns elements and networks in 
which the followingtwo rules apply: (Kirchoff's laws) 


1. Thealgebraic sum of the node values around 
any closed path (loop) in a network is zero. 


2. Thealgebraic sum of the transmission quan- 
tities at any node is zero. 


: Korn and Huskey, "Computer Handbook", Chapter 2, 
McGraw-Hill, New York, 1962. 


TABLE I. Single-element Computing Blocks 


PHYSICAL SCHEMATIC MATHEMATICAL 
ELEMENT SYMBOL RELATIONSHIP PROGRAM BLOCK 


Multiplier 
(POT) 
Y = KX 


Inverter 
= 
Y=-X 


Integrator 
y-> 
Ss 


Differentiator 
Y=sxX 


Resistor 


Torsional 
Spring 


Heat 
Conductor 


Viscous 
Damper 


TABLE I (Continued) 


PHYSICAL SCHEMATIC MATHEMATICAL 
ELEMENT SYMBOL RELATIONSHIP PROGRAM BLOCK 


Electrical 
Capacitor 


Viscous 
Damper 


Mechanical 
Mass 


Thermal 
Mass 


Electrical 
Inductor 


Mechanical 
Mass 
COMPUTER PROGRAM MECHANIZATION a) assigning aposition inthe mechanization 
COMMENTS diagram to all necessary variables and 
operators. 


b) accomplishing the necessary intercon- 


1. Operational amplifiers operate with the junc- nections so that each output has its pro- 
tion at virtually zero potential; the currents per inputs. 
into the junction through all the input and 
feedback paths must sum to zero; and there c) checkingto see that all input/output rela- 
is asign inversion between the amplifier in- tionships in the mechanization are ‘satis- 
put and output quantities. fied. 
2. Anyclosed-loop ina problem mechanization EXAMPLES OF NETWORK PROGRAMMING 
Simulating a stable physical system will gen- 
erally have an odd number of sign inversions. The following examples illustrate some of the tech- 
niques used in applying the block programming meth- 
3. In mechanizing an equation, set the highest ods. Examples are given for linear elements only. 
order derivative only on the left of the equals Non-linearities can be included simply by replacing 
sign to obviate the necessity for differen- the potentiometer corresponding tothe non-linear ele- 
tiation. ment with a multiplier, function generator, or other 
suitable computing element. Backlash, hysteresis, 
4. The mechanization is accomplished by: stiction, etc can be included as required. 


TABLE Il. Two-Element Computing Blocks 


MATHEMATICAL 
SCHEMATIC RELATIONSHIP PROGRAM BLOCK 


MASS - SPRING - DAMPER SYSTEM 


XV Xo1Vo 


MOTION IN X 


F DIRECTION ONLY 


Element - force equations 


Spring K,;: F., = K, | 1 


Spring Ky: Ps -XK Initial values att = 0. 


X, 0), X_ 0), v, 0), vy (0) 


) 


Damper B: Fy = B (Vo - Vv 
Forces acting on masses 


F ~ Fey r F 59 D 


m2 ~ ~“Fs2 ~ Fp + F 


Steps to mechanize 


1. Assume X,; -X,, “Vy> Vo and F are available 


2. Draw blocks for M)> My 


3. Draw blocks for K,, Ky, B 


4. Make input connections to mass blocks to produce Xi, -Xo, -V 


———-— — —  - 
Pema | tt pe 


| 
| 
| 
| 
| 
|" 
| 
| 
| 
| 
| 


SPRING! | MASSI| |DAMPER ! SPRING 2 MASS 2 


ESR ARE Rea ees eerie ees Lemur meres 


The above mechanization makes each parameter avail- 
able as an individual adjustment. The programming 
is easily accomplished by generating the forces im- 
posed on the elements and then operating on the forces 
to obtain velocity and displacement. Any or all of the 
initial conditions, X,(0), X2(0), vz(0) and vg(0), may 
be zero, or may also be variable. 


Note that the block representing the damper contains 


two summers. The first summer produces (vy - vg), 
and the second summer acts only as a signchanger. 
The two-amplifier, one-pot combination can be reduced 
to a more simple two-pot combination by connecting 
V1 and vo through B pots to the mass block inputs. 
The less complex alternate isshown below. This con- 
figuration requires fewer amplifiers, but variations 
in M,> Mo, or B necessitate two or three adjustments. 


* Do not connect to “Vy to get pot setting of B instead 


of =. A coefficient pot output cannot be connected 
1 


to another coefficient pot input. This is a practical 
limitation of the computer hardware. 


Rigid Bar Supported on Springs 


e 22 11 F 
§= 
sI 
F, = -K,Y, 
Fo = -Ky Yo 


Rigid bar of mass M, momentof inertia about C.G. of 
I. Assume 6smallsuchthat 6=sin6. Vertical motion 
only. Yo is displacement of C.G. 


a F,+F,+F 
0-7 sM 

Y¥, = Yo- £48 

Yo = Yo + £99 


MASS 


SPRING 2 [ROTATIONAL | 


_ 
| 
| 
aE 


SPRING | 
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| 
aren eee 


One-Dimensional Heat Transfer 
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M-MASS 


C-SPECIFIC HEAT 


Electrical R-L-C Networks 


RESISTOR | aeileectiepcsl Wace CAPACITOR Mesa RESISTOR | 


Sane ee ec) RR Eee Lee nney ayeecneny es 


Io Iz 
1, I, 
ee eae Ba Bs = 3C, BE, = 1, Rg 
E 
foe oe ee 


TWIN-T NETWORK 


1 3 3 2 
I. = -E, = 
2 R, 3 SC 
ee ar at 
3 Ro 
-(I, +I,) 
8 3° 
(E, - E,) = sC, 
I, - I, 
(Ey 7 E,) ~~ §C 


E, = (I, - I, +1, + Ig) 


E, = (E, - Ey) +E; E, = Ey - (Ey - E,) 


-(T,-Io+I3+Ig) 


C)R3 
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SCALING 


The mechanization diagrams associated with the pre- 
ceding examples show the computing element inter- 
connections without indicating potentiometer settings, 
amplifier input resistor values, andintegrator capac- 
itor values. These are unscaled mechanizations. 
Scaling is required to relate the computer voltages to 
the physical problem variables. 


Tne analog computer is voltage limited to a practical 
range of operation of +100 volts as a maximum and +1 
volt or zero aS a minimum. It is usually not possible 
tonumerically equate analog voltages to physical vari- 
ables. That is, the magnitude of physical variables 
must be scaled to fall within the useful range of the 
machine. This procedure is called amplitude scaling. 


The analog computer is also limited as tothe speed with 
which it will solve a problem. Practical bounds on 
solution time are 100 seconds as a maximum and 100 
milliseconds as a minimum. Events in the physical 
world usually occur in time intervals which fall out- 
side these limits. Thus, the simulation of real world 
phenomena with the analog computer is ordinarily 
faster or slower. Theprocedure for relating computer 
time to physical time is called time scaling. 


Correct scaling is an important factor in reducing sim- 
ulation errors. An error analysis of any particular 
program is about as complex as the problem which is 
being solved by the program. Consequently, except 
for unusual cases, an analysis of error is not made. 
An estimate of the accuracy of the program can be 
found from check cases for whichthe answer is already 
known. Error is reduced by scaling the program so 
that all potentiometer values and amplifier gains are 
reasonable. Frequently, in the process of scaling a 
problem, it will be found that some parts of the model 
are not significant and can be eliminated. Scaling is 
a good check on reasonability. 


Scaling is anart, not a science. A good deal of exper- 
ience is required to become proficient. Scaling canbe 
done inseveral ways. Thetechniques discussed in this 
chapter are those in most prevalent use among analog 
programmers. 


If X is a physical variable corresponding to computer 
voltage, V, then ascale factor, a, is chosen so that 


V=aa-X 


will fall within the practical range of the computer. It 
is good practice to make V as large as possible. If 

is physical (real) time corresponding tot,, the com- 
puting time, then these two can be related by a scale 
factor, N: 


t = Nt, a” = N™a™ 
Cc js a p 


It follows that 
S(F)at,, = Nf (FUMt,,. (1) 


N is chosen so that the problem solution time on the 
analog computer falls within reasonable limits. Prac- 
tically, this is accomplished by a choice of N which 
will not result in prohibitive amplifier gains. 


The general procedure for scaling is: 


1. Generate anunscaled program having a math- 
ematical structure which agrees with the pro- 
blem (model) to be solved. 


2. Identify the location of all voltages corres- 
ponding to physical variables. 


3. Associate a scale factor with each of these 
voltages (e.g. vo = aX). 


4, Label the program with the scaled physical 
variables (i.e. aX) rather than the voltages. 


5. Estimate maximum values of the physical 
variables. 


6. Choose the amplitude scale factors andthe 
time scale factor, N, so that all pot settings 
andgains are reasonable. 

In order to accomplish step (4) it is necessary to know 
how input and output scale factors are related for 
various computer elements. 

POTENTIOMETER 


Physical equation: Y = bX 


=a X, V_=a_Y 
Scale factors Ve y y 
a 
Scaled equation Vy = ae Ve 
Program aX ( ) he 
ba 
a 
a 
x 


Note that ascale change can be made withapot. Assum- 
ing either a, or ay has been established previously, 
the other is choseh so that ba,/ay is a reasonable pot 
setting. The potcan be used solely fora scale change. 
In that case b = 1. 


SUMMER 
Physical 
equation: Z = X + Y 
Scale 
factors a5 = a2, we = aX, Vy = 2° 
Scaled a, a, 
equation: vy “ai Vy + a vy 
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Now, if only one-gain inputs for the summer are used, 
then it must be true that a =a, =a and the pro- 
gram is y 


-a,X 


a7Z 


Thus, it is seen that input scale factors must be the 
same for a one-gain summer and that a change of 
scale cannot be made through asummer. However, 
a different input scale factor can be used fora10-gain 
input. The program for this case is, 


-a,X 


< a,Z 
ayY 0 


where it must be true that a. = 10a... 


INTEGRATOR 


As with the summer, the input scale factors (except 
for the 10-gain input) must be the same. However, a 
change of scale can be made with an integrator. Both 
amplitude and time scaling can be done. Scaling for 
an integrator is derived below, where C is the value 
of the capacitor in ywfd, and R is in megohms. 


Physical 
equation: Y =f Xat 


Scale 

factors Vi= a_.Y, V_=aX,t._= Nt 
y y x <x Cc 

Scaled ay 1. 

equation: Vy = a, W Vit. 


Thus, the integrator gain, ao must be 


RC 


ts dat ad 
~ ay (2) 


and the program is 


— A,X NayY 
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MULTIPLIER 
Physical 
equation: XY = 2 
Scale 
factors: Vo = aX, V_=a Y, Vi_=aZ 
x x y y Z Z 
Scaled 
equation: V_V 100a 
xX y Zooy 
100 aa Z 
x y 
Program: 
OxX a,Z 
ayy 
It must be true that 
100a 
Z 
= 1. 
aa 
x y 


From this requirement itis seen that the relationship 
among the scale factors is the same as thatfor input and 
output multiplier voltages: 


a = <~V 
“100 
DIVIDER 
Physical X 
equation: = = Z 
a 4 
Scale 
factors ver aX, vy = a¥, Ne = a, Z 
Scaled Vs aa. 
equation: 1005— : 00. = NS 
y x 
Scale 
factor a, 
constraint: BF 100-— 
y 
Program: 
a,X a,Z 


SQUARE ROOT 


Physical 

equation: Y = ~X 

Scale 

factors v= ayY, V5 = ax 

Scaled a 

equation: V, = y . /V, 
10 /a,, 

Scale 

factor 


constraint: a, = 10 fa, 


Program: 


~D-» 


FUNCTION GENERATOR 


Physical 
equation: Y = f(x) 


Scale 
factors V = ayy, ve 


Scaled 


equation: V_= a f(—~* 
y y a 


EXAMPLE: 


Mass-spring-damper system 


dx 
ome 7X 
dtp 
Physical 
equations: Spring Force, F_ = 


Ss 


Damper Force, F, = 


D 


Mass 2 -(F_+F,) 
Acceleration, o* Bei ees 
dt 

p 
, -4 
M=1;B=4x10 


dx 


—. (0) = 0, X (0) = .95 


dt 
p 


>K=6.4x10 


dt 


Estimates: (ee) =8x 107°; Xmax. = 1 
p 


The unscaled program is: 


The scaled program is obtained by assigning literal 
scale factors, determining their numerical values, and 
calculating the coefficient potentiometer settings. 


4 


p------- 


1 


The -> and + factors appear 
N2 N 

due to the inherent gain of N 

through each integrator, from 


Eq. (1). 


In general, the determination of the numerical values 
involves atrial and error approach. Any convenient 
parameter can be picked for the starting point. Ifa 
selected value leads to unreasonable potentiometer 
settings or unreasonable amplifier gains, then new 
values may be required. Reasonable pot gains are 
between .1 and1. Reasonable amplifier gains are 
between 1 and 20. Amplifier gains should be integer 
values. Non-integer values can be obtained from a 
potentiometer, in combination with gain-of-10 inputs, 
if necessary. 


1. Determination of a. In order to utilize the 
full dynamic range of integrator 2, set 
aX max = 100 (volts). 


a = 100 
13 
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Potentiometer 1 setting. 


a X (0) 


Pil = 


Determination of a At this time, a trial 


value for N maybe selected. Analog solution 
times of .1 to 100 seconds, or computer nat- 
ural frequencies from .1 to 100 radians per 
second are the preferable operating ranges. 
Anestimate of the physical time or frequency 
ranges can be used to select a value for N. 


For this example, a valueof N = .01 appears 
reasonable. 


b _ 100 (volts) 


N ~ (dX/dt, max, (4) 


Looking ahead, note that integrator 2 will 
aN 
— 
herent in the integrator from Eq. 1. There- 


require a gain of The gain of N is in- 


fore, the = value for the integrator will be 


= It is highly desirable that this factor be 


aninteger, preferably either 1 or 10. There- 
fore, b = 100 will be used. The effect of this 
is to force integrator 1 to utilize something 
less than its fulldynamic range, but the effect 
is not serious inthis example. As a general 
rule, it is not possible to scale so that all 
amplifiers work over their full range. 


4. 


Potentiometer 2 setting; 


P2 = x = (6.4x 107 %)c (5) 
N aM 


The value of c ischosento assure that a rea- 
sonable value for P2 is obtained and to assure 


that — is an integer. 
c = 100 will satisfy both conditions. 


P2 = .64 


Potentiometer 3 setting: 


cB _ 04 


P3 = NoM © 


This setting for potentiometer 3 is below the 
desired minimum potentiometer setting value 
of .1. Inthis example, however, it is not pos- 
sible to have a large setting for potentiometer 3. 
This follows from a consideration of the phy- 
sical problem. The damper, B, was chosen 
tohavea very slight effect on the mass-spring 
system, so it necessarily follows that the 
portion of the analog program related to the 
damper willhave onlya slight effect. It should 
be noted that abnormally low potentiometer 
settings will frequently arise in simulations of 
systems containing negligibly small elements. 
It should be recognized that the low settings 
stem from the actual system characteristics. 
The programmer, therefore, need not spend 
time in futile attempts to improve the scaling. 


The final program with numerical values is: 


ee 1074 Vi 
dtp 


2 
d ».4 = 
420. \Yp2*¥ps) 
p 


GENERAL COMMENTS ON TIME SCALING: 


Tne following points concerning time scaling are of 
particular importance: 


A, 


If all integrator gains are simultaneously 
changed bya factor k, the time scale changes 
from N to, and all computer frequencies 
increase by k. The amplitude scaling, how- 
ever, does not change. All numerical scale 
factors, all potentiometer settings, and all 
computer problem voltages remain the same. 
This property of amplitude invariance under 
time scale change allows the programmer to 
easily speed up or slow down the computer 
solution. 


The most convenient way of changing all inte- 
grator gains simultaneously is by changing 
the capacitors by a fixed factor. Integrator 
capacitors can be simultaneously changed by 
means of the computer Time Scale control-. 


An event which occurs in a fixed real-time 
interval can occur in several different com- 
puter time intervals depending on the choice 
of integrator capacitors, that is, depending 
on the time scale factor. 


Events of identical character which occur in 
different real-time intervals can all be made 
to occur in the same computer time interval 
by appropriate time scaling. 
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